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Atomic Bose-Einstein condensate is heated by atomic losses. Predicted depletion ranges from 1% 
for a uniform 3D condensate to around 10% for a quasi-lD condensate in a harmonic trap. 
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Heating by atomic losses. — An ideal Bose-Einstein 
condensate (BEC) is a state where all bosons occupy the 
same single particle state J>q. So far the atomic BEC 
is the closest to this ideal [if. However, even the atomic 
BEC is not perfect because atoms are depleted from 4> by 
atom-atom interactions (quantum depletion) and ther- 
mal fluctuations (thermal depletion). Moreover, atomic 
condensates last only for tens of seconds before atomic 
losses empty the trap. This process is dominated by 3- 
body losses, where 3 atoms collide to form a bi-atomic 
molecule and an atom with large kinetic energy, and then 
both of them escape from the trap. There are also 1-body 
losses, where individual atoms are kicked out from the 
trap by external heating agents, see e.g. Ref.|2|. While 
the 1-body losses can, at least in principle, be minimized, 
so far the more intrinsic 3-body losses have not been elim- 
inated. 

Atomic losses from a trap can be modeled by repeated 
action of an annihilation operator ip(x). An ideal con- 
densate \N : (f>o) is a "fixed point" of the annihilation 
operator: action of V>(x) neither depletes an ideal con- 
densate nor changes its wave function </> , ip(x)\N : <fi ) — 
\fN4>o(x)\N — 1 : 0o )■ A destructive measurement of 
atomic positions in a trap can also be described by re- 
peated application of $(x). In Refs.|3] it was shown on a 
few examples that such a repeated annihilation can grad- 
ually increase a condensed fraction of atoms remaining 
in a trap. Annihilation is driving the remaining atoms 
toward an ideal condensate. This phenomenon is a foun- 
dation for a quite popular believe that "atomic losses 
improve a condensate" . This is not a sound foundation. 

In the destructive measurement all annihilations hap- 
pen at the same time, or in a very short measurement 
time. Hamiltonian of atoms has no time to do anything 
before all atoms are annihilated. In contrast, in a typical 
experiment it is the Hamiltonian that is much faster than 
atomic losses. An ideal condensate \N : 4>q) is not an 
eigenstate of the Hamiltonian: atom-atom interactions 
are depleting atoms from the condensate wave function 
<f>o ■ On the other hand, eigenstates of the Hamiltonian H 
are not "fixed points" of the annihilation operator: ip(x) 
applied to an eigenstate of H gives a non-stationary state 
because ^(x) does not commute with H, [^j(x),H] ^ 0. 
A competition between the Hamiltonian depleting a con- 
densate and the atomic losses increasing the condensed 
fraction leads to a state that is neither an ideal con- 



densate nor the iV-body ground state. A condensate is 
heated by atomic losses — this phenomenon will be more 
quantified in the following. 

In Ref.0 Timmermans demonstrated that heating of 
fermionic atoms by atomic losses is a serious obstacle 
on the way to atomic superconductors. In the case of 
fermions the mechanism is simple: atomic losses create 
holes in the Fermi sea. In the case of bosons the heating 
is a conceptually more subtle effect. One has to realize 
the interplay between the losses and the Hamiltonian to 
see that losses are in fact heating a condensate. However, 
as we will see below, a conceptually subtle effect is not 
necessarily a quantitatively subtle effect. It is a serious 
and rather fundamental limitation on quantum coherence 
of atomic BEC. 

Master equation. — The non-unitary (due to atomic 
losses) evolution of trapped atoms is described by a mas- 
ter equation || 
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Here H is a Hamiltonian of trapped atoms 



H = I d A x ^-V^VtA + V(x)$4> + 1^7/^ 
2m 2 

' (2) 

with g — AttTi a/m, where a is the s-wave scattering 
length. V ip l (x) is a Lindblad superoperator |fj, 



V[a]p = a pa* 



1^ f 
-a'ap- -pa*a 
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describing ^-body losses. 

Bogoliubov theory. — We assume that almost all N 
atoms occupy a common condensate wave function 4>o (x) , 
which solves a stationary Gross-Pitaievskii equation |jj 
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The annihilation operator can be split into a condensed 
part and a non-condensed part which is then approxi- 
mated by an expansion in Bogoliubov modes 0, 



ip(x) 



OO 

A^ 1/2 O (x)+ [brnU m (x) 
m— 1 



y m {x) .(5) 
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Here b s are bosonic quasiparticle annihilation operators, 
and the wave functions u rn and v m satisfy Bogoliubov-de 
Gennes equations ||- 

The operator of a number of atoms depleted from <p 
is dN = J d 3 x dip'' 'Sip. In the Bogoliubov vacuum state 
1 0&) without any quasiparticles, 6 m |0b) = 0, the number 
of depleted atoms is 



^ (0) = E j d " x m 2 ^E^™- 



(6) 



More generally, in a state with exactly n m quasiparticles 
in a mode to the number of depleted atoms is 



dN = J2 [dN™ + (l + 2dN^) n ri 



(7) 



Here we used J d 3 x (|u m | 2 — l^ml 2 ) = 1. 

Master equation in the quasiparticle 
representation. — Bogoliubov expansion (J5J can be 
used to rewrite the master equation Q in the quasipar- 
ticle representation. Expansion to second order in 6's in 
Eqs.J2CJ results in a Bogoliubov Hamiltonian which 
is a sum of harmonic oscillators H w Yl m ^rn&m^mj 
and an approximate master equation 



dp 



fa = E -iUrnPrJm, P\ + 



(8) 
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(1 + ni m )V[b m ]p + n lm Vp m ]p 

The coefficients are defined by integrals 

d 3 x \<j> \ 2l - 2 \u m \ 2 = a lm (l + n lm ), (9) 

d 3 x \<j) \ 2l ~ 2 \v m \ 2 = ai m ni m . (10) 

In addition to small depletion, derivation of the master 
equation (jSJ) requires the rotating wave approximation 
(RWA). In the RWA we neglect all terms of the form bb, 
or b m b\ for to ^ n, but keep all terms like &m6 m or 
b] n b m . The RWA is accurate when the Hamiltonian evo- 
lution is much faster than atomic losses, or more precisely 
w m 3> 'Yliil"fiaim,N l ~ 1 ni m . This condition is satisfied in 
all present day experiments. 

A thermal state. — Due to atomic losses the coeffi- 
cients in IJSJ are not constant. However, we fix them (for 
a while) and analyze a stationary state of the resulting 
master equation. Later on we will see that such an anal- 
ysis allows us to predict a lower bound for a stationary 
depletion of a condensate caused by atomic losses. 

A remarkable thing is that the master equation JSJ, 
with the coefficients fixed, can be recognized to describe 
a set of harmonic oscillators (numbered by to) coupled 
to external heat reservoirs (numbered by I). Every oscil- 
lator relaxes to a thermal state. When atomic losses are 



dominated by only one of the channels I, then average 
numbers of quasiparticles in the thermal states are 



= Tr tyj>m P(t ~> 00 ) = 



"I, 



(11) 



When many channels I are involved, then the averages 
n m can be obtained from equations 
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T,ih«*imN l - 1 (l + ni m ) 
J2i hiai m N l - 1 ni m 



(12) 



Every oscillator to can be assigned to an inverse temper- 
ature P m which follows from a textbook formula n m — 
( e f3 m hu m _ Thc thermal state is 



p(t 



oo 



(13) 
For a 3D 



dN™ are 



A thermal state of a uniform BEC 

condensate with <po(x) — const, the n; m 
independent of I, compare Eas. (|9llOI6[l . In a uniform 
condensate of density p c a phonon of momentum hk has 
energy 0] 



Hulk = 



in 2 k 2 (h 2 k 2 
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Here c = y/ ' gp c /m is a velocity of sound. In the thermal 
state (fT~5f) there are on average 



n k = / d 3 x \v k \ 2 = 
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phonons of momentum hk. With the general formula 
we can calculate a fraction of depleted atoms in the 
thermal state d = dN/N — yj Amp c a 3 . For a typical 
condensate density of p c — 10 20 m~ 3 we find depletions 
^ 23 Na = 0.44% and c?87 Rb = 1.55% at the scattering 
lengths of a = 2.5 nm and a = 5.8 nm respectively. 

The depletion is dominated by a contribution from 
small k where the number of quasiparticles n k in Eq. i|15|) 
is divergent. A remarkable thing is that for small k the 
equipartition of energy n k huj}~ — (3~ l yields the same 
temperature T for all phonons. The temperature is 
T23 Na = 76 nK and T«7 Rb = 37 nK. 

These estimates are valid for a 3D uniform conden- 
sate. In less than 3D the infrared divergence n k ~ k^ 1 
results in a divergent depleted fraction d. Anticipating a 
much larger but finite d we now turn to effectively one- 
dimensional harmonic traps. 

A thermal state of a BEC in a ID harmonic 
trap. — In a sufficiently anisotropic trap V(x,y,z) = 
^muj 2 [x 2 + n 2 (y 2 + z 2 )] with n » 1 the y — z state of 
all atoms is frozen in the ground state. The condensate 
wave function 4>o(x) solves a ID Gross-Pitaievskii equa- 
tion (@J with an effective gm — g/S, 2 , where £ is the size 
of the ground state in the y — z plane. We solved the 
ID Bogoliubov-de Gennes equations @ to get u m and 
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TABLE I: Lower bound for the stationary depletion due to 
atomic losses for parameters corresponding to quasi-lD |t| 
and ID [Toil experiments. In the calculations the values of 
gioN/Tiui have been estimated to be 7500 (where N — 1.5 ■ 
10 5 ) and 500 (N = 10 4 ), respectively. 





1-body losses 


3-body losses 


quasi- ID 


6% 


10% 


ID 


2% 


4% 



Vjn for two sets of par ameters relevant to the quasi- ID 
[9j and strictly ID [lOj experiments. Relative depletions 
d = dN/N corresponding to the thermal state (|13|) are 
listed in Table [I] These values give lower bounds for the 
stationary depletions in these experiments. The 10% de- 
pletion in the 3-body losses for the quasi-lD condensate 
is actually close to the thermal cloud fraction estimated 
in the Hannover experiment 0- In other words, this ex- 
periment is close to the minimal stationary depletion set 
by atomic losses. 

Relaxation time. — It takes time to reach the ther- 
mal state l|13f) . To simplify notation we assume here 
that losses are dominated by only one of the channels 
I. Time evolution of the average number of quasiparti- 
cles n m (t) — Trp^b^bm directly follows from the master 
equation JHJ) and satisfies a differential equation 



dt 



n m (t) = -IjiaimN (t)[n m (t) - m 



(16) 



n m (t) is relaxing toward its equilibrium value ni m . 
Ea. (|16|l has to be compared with the decay law for the 
total number of atoms N(t) = Tiyj(i) J d 3 x ifi^ip 



dt 



N(t) = -Z 7 ,a z W'(t) . 



(17) 



Here on = J d 3 x \4>o\ 21 - This equation is valid for small 
depletion, when almost all atoms occupy 4>q(x). 

As dominant Bogoliubov modes are localized on the 
condensate, we can approximate ev; m » a; . Consequently 
the relaxation and decay rates in Eas. (|16|) and (|17fl are 
comparable, dn m /n m w dN/N . As the equilibrium value 
m m in Ea. l|16|l depends on the time-dependent N(t) and 
the rates are comparable, n m (t) will never quite reach the 
instantaneous equilibrium value ni m (t). In the following 
we will consider 1- and 3-body losses in the two mode 
(double well) toy model to see that the estimated equi- 
librium depletion, we have considered so far, is a lower 
bound for a stationary depletion. 

Double well model. — The model is described by 
Hubbard Hamiltonian 



H2 = -fl{d\d2 + aldi) + -(d\a\didi + a\a\a-2.a2) 



(18) 



The master equation Q becomes 
dp 1 



dt 



-[H 2 ,p] 



4=1,2 



(19) 



We assume that Z-body losses dominate. 

The condensate wavefunction (f> = (1,1)/ y/2 and 4>\ — 
(1, — l)/y2 span the single particle Hilbert space. There 

= (1 



is only one Bogoliubov mode with u± 
and v\ = —6\I\JX 2 — 1, where X 



4n \ 

N . 



(l + — 1- As </>o is "uniform", the equilibrium 
number of quasiparticles is independent of I, 



n n = dN[ a) = v\ Vl = ^ \ + 0(N-^) , (20) 
Depletion in the equilibrium thermal state (|13J) is 



(0) 



1 + 2dN{ 



(0) 



N 



1 
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+ 0(N~ 1 ) . (21) 



1-body losses. — We begin with I = 1 and a conden- 
sate initially in the Bogoliubov vacuum state, i.e. the 
number of quasiparticles ni(0) = 0. A formal solution of 
Eci. (|lt)|) . where an = 1, is n\{t) = 71 J Q dt' exp[— 71 (t — 
tf)]nn(t'). N(t) = N e~ 7lt , which solves Ea.^T7)l. and 
Eas. H2QI7|) give (time-dependent) depletion 



d l=1 (t) = 



W7 



f 



V327T/V 



O^N- 1 ) 



(22) 



Here / = N(t)/Na is a fraction of atoms remaining in 
the trap. When / — * (the system is forgetting about 
the initial conditions) the depletion ll'L'l) becomes roughly 
twice the equilibrium value l|21|). At small / the system 
is reaching a stationary state with twice the equilibrium 
depletion. 

When we initially prepare the system with ni(0) = 
2\J — |, instead of the rather arbitrary ni(0) = 0, 
then the depletion is stationary — it does not depend on 
/ (and through / on the initial No), 



•vJ=i 



1 



oiN- 1 ) 



(23) 



This stationary depletion 112MI) is roughly twice the equi- 
librium value l|21|) . 

3-body losses. — When 3-body losses dominate, the 
system gets much closer to the equilibrium than in the 
case of 1-body losses. For a condensate initially in the 
Bogoliubov vacuum state, ni(0) = 0, a similar procedure 
as for / = 1 leads to a depletion 



jZ=3 



(*) 



3(1 - f 5 / 2 ) 



+ 0(N- 1 ). (24) 



40O V32QN 
When / — > d l=3 (t) is approaching a stationary value 



Here we use rescaled dimensionless units such that h — 1 . 



(25) 
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that is only 1.2 higher than the equilibrium value (|21(l . 
Starting with the initial ni(0) = |^/ 5^ ~ \ resu lts in 
a depletion independent of / (and thus also on No) and 
equal to the stationary value l(2*B|) . 

We conclude that our estimates of depletion based on 
the equilibrium values are lower bounds for stationary 
depletions. 

Numerical experiment. — We verified the predic- 
tions I|22I23I24I25|) of the Bogoliubov theory in numeri- 
cal simulations. For a large No a direct solution of the 
master equation l|19|) is not the most efficient. It is bet- 
ter to replace the deterministic p(t) by an ensemble of 
stochastic pure states \^(t)), such that p(t) is repro- 
duced as an average over many stochastic realizations, 
= p(t). A stochastic "unraveling" of the 
master equation l|19|) is given by Ito stochastic nonlinear 
Schrodinger equation [llj 



dm 



-idtH 9 m 



dtji 
2 



E 

j=l,2 



E ^(*) 



.7=1,2 



*5 



I*) 



(26) 



where (A) = <Wj(t) £ {0,1} is a stochastic 

process [dAfj(t) = 1 when I atoms escape from a well j 
between t and t+dt, and otherwise]. A probability that 
I atoms will escape between t and t + dt is df((a] )'q'). 
In FigQ]we compare the predictions (|22l23l24l25|li with 
corresponding averages over many stochastic realizations. 
As they compare quite well, the Bogoliubov theory © 
passes the test on the double-well model. 

Conclusion. — A condensate is heated by atomic 
losses. The depletion of the system approaches a sta- 
tionary value that ranges from around 1% for a uniform 
3D condensate to around 10% for a quasi-lD harmonic 
trap. As atomic losses cannot be easily eliminated, this 
depletion is a serious limitation on quantum coherence of 
atomic BEC. We only note here that outcoupling in the 
atom laser is a non-markovian process of atomic losses. 
Its influence on laser coherence will be addressed else- 
where. 
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FIG. 1: Relative depletion d = dN/N in the 2-mode model, 
as a function of the number of atoms that escaped from 
the trap due to 3-body losses (for No = 500, Q. = 0.5 and 
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results). The dotted line is the equilibrium value (12 U . The 
inset shows the same as in the main figure but for 1-body 
losses with 71 = 10~ 2 . 
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